Temperature dependent spin susceptibility in a two-dimensional metal. 
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We consider a two-dimensional electron system with Coulomb interaction between particles at a finite tem- 
perature T . We show that the dynamic Kohn anomaly in the response function at 2k F leads to a non-analytic, 
linear-in-J correction to the spin susceptibility, 6x(T) = AT, same as in systems with short-range interaction. 
We show that the singularity of the Coulomb interaction at q = does not invalidate the expansion of A in pow- 
ers of r s , but makes the expansion non-analytic. We argue that the linear temperature dependence is consistent 
with the general structure of Landau theory and can be viewed as originating from the non-analytic component 
of the Landau function near the Fermi surface. 



Introduction — There has been substantial recent interest in 
the temperature dependence of various Fermi liquids proper- 
ties for both short-range and long-range interactions between 
particles The revival of interest in the prob- 

lem is two-fold. On the experimental side, technical advances 
now allow one to measure the temperature dependence of the 
thermodynamic parameters such as specific heat and spin sus- 
ceptibility in "classical" 2D Fermi liquids with short-range 
interaction, such as monolayers of 3 He, as well as study two- 
dimensional semiconductor structures with long-range inter- 
action and with relatively low Fermi temperatures (~ 1 K). On 
the theory side, the leading interaction corrections turn out to 
be non-analytic functions of temperature making the subject 
particularly important. 

Nai've power counting arguments suggest that the temper- 
ature dependence of any thermodynamic quantity, includ- 
ing the spin susceptibility and the specific heat coefficient 
C(T)/T = y, should start with terms quadratic in temperature. 
This conjecture is based on the observation that a thermody- 
namic quantity at a finite temperature typically can be written 
as J a(e)n(s)ds, where n(e) is the Fermi distribution function 
and a(e) is some function. If the latter is smooth, the temper- 
ature dependence starts with a term of order T 2 |7]. Such a 
temperature correction is called "analytic." This is also con- 
sistent with the intuitive expectation of the one-to-one corre- 
spondence between the non-interacting Fermi gas and the in- 
teracting Fermi liquid since in the Fermi gas, the Sommerfeld 
expansion leads to simple quadratic temperature corrections. 

However, the assumption about the analyticity of the func- 
tions involved in the calculation of various thermodynamic 
properties of the Fermi liquid is quite generally not justified 
because in any Fermi liquid, the dynamic interaction between 
particles gives rise to a non-analytic energy dependence of 
a(e). This leads to temperature corrections which do not scale 
as T 2 and are therefore called "non-analytic." Collecting these 
non-analytic corrections is a subtle theoretical problem. 

The subject of this paper is the temperature correction to the 
spin susceptibility for a 2D system of fermions interacting via 
a long-range Coulomb interaction. For Fermi systems with 
short-range interaction, perturbative calculations for a model 
with a small U(q) have demonstrated that 5xAT) is linear in 



T and that the prefactor depends only on U(2pf). This re- 
sult, however, was obtained under the assumption that for all 
q, the dimensionless Born parameter u(q) = mU(q)l(2n) is 
small. For Coulomb interaction this is obviously not the case, 
and one has to verify explicitly whether the linear dependence 
8Xs(T) °c T still holds, and whether the prefactor can be ex- 
panded in r s . The spin susceptibility has recently been mea- 
sured at various T in Si inversion layers iill9[ flolfTlll . and a 
quantitative theory is required to interpret the temperature de- 
pendence of the experimentally measured^,. In what follows, 
we compute the spin susceptibility in the perturbation theory 
and beyond, and relate the prefactor of the linear in T term in 
an arbitrary Fermi liquid to the spin component of the quasi- 
particle scattering amplitude at the scattering angle 9 — n. 

We also consider in detail the relation between the non- 
analytic T dependence of the thermodynamic parameters and 
Landau Fermi liquid theory. The Landau theory operates with 
the quasiparticle interaction function for the particles at the 
Fermi surface. In this theory, the spin susceptibility is inde- 
pendent of T, and is expressed via the particular partial com- 
ponent of the Landau function. The temperature corrections 
to the Fermi liquid theory come from quasiparticles which 
are slightly off the Fermi surface. We demonstrate explic- 
itly that non-analytic corrections to the spin susceptibility can 
be viewed as originating from the non-analytic momentum 
dependence of the quasiparticle interaction function f(p,p') 
at small deviations from the Fermi surface. We argue that 
f(p,p') is non-analytic in deviations from the Fermi surface, 
and that this non-analyticity gives rise to the emergence of the 
linear in T terms in the gr-factor and the spin susceptibility. 

Perturbation theory — Consider first the case when the 
Coulomb interaction is weak at q ~ pf, i.e., when r s <k 1. 
The temperature correction to the spin susceptibility can be 
calculated either by explicitly evaluating the static particle- 
hole polarization bubble in a zero field, with insertions due 
to the interaction, or by evaluating the free energy in a finite 
magnetic field H and then differentiating over H. Either way, 
one obtains that, to the leading (second) order in the interac- 
tion, the linear in T term in the spin susceptibility comes from 
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2p p processes, and 

SXs(T) = ^p au ii (^) 2 |-, for T/E F « 1, (1) 

where ^p au ii = fyn/n is the susceptibility of free fermions, 
and r s - ^j2mU(2p F )ln. The special role of2pf terms in the 
perturbation theory can be most easily understood by eval- 
uating 6xs(T) from the free energy. To second order in the 
interaction, the free energy Q2 consists of two particle-hole 
bubbles connected by U(q) 



(2) 
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where a, {3, y and 8 are spin components of four fermions in- 
volved. Since the Coulomb interaction is spin independent, 
the spins of the two fermions within each bubble are parallel 
(i.e., a — B and y = 6). However the spins in different bubbles 
can be either parallel or antiparallel to each other. The singular 
6xs(T) comes from the antiparallel spin configuration between 
bubbles. For such spin orientation, fermions in each bubble 
have different Fermi momenta p F and p F due to the Zeeman 
splitting. Near 2pf, this splitting is relevant as the polariza- 
tion bubble is non-analytic in both momentum and frequency. 
The non-analytic momentum dependence is normally associ- 
ated with the Kohn anomaly and related Friedel oscillations. 
For the T dependence of the spin susceptibility, however, one 
actually needs the dynamic polarization bubble. In 2D, the 
singular part of the polarization bubble behaves near 2pf as 



i lu , o s .: f { J. + i 9s\ 2 - i + 0__^) 2 _ L(3) 



At small frequencies and q < 2p F , this reduces to 



FL(q, Q m ) oc 
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y/2p~ f 



(4) 



Integrating the product Yl ++ (q, O m )II (q, D. m ) in (0 over q, 
we find that the frequency dependence is not analytic: Q2 K 
U 2 (2p F )T En„, O^ log [Q 2 + (fi B H) 2 ]. Evaluating the sum 
one finds that O2 contains a cross term TH 2 II 311 . Differen- 
tiating over frequency, one then obtains 5x s (T) « T, as in Q. 

At the same time, the potentially dangerous small q region, 
where the Coulomb interaction is large, does not contribute 
to Q (in this respect, our results differ from those in Ref 
0). The reason is that at small q each of the two polariza- 
tion bubbles contains only a non-singular, multiplicative de- 
pendence on the magnetic field; in two dimensions, this de- 
pendence comes through v F (H) - pf(H)/m in Tl(q,D. m ) - 

(m/2n) 1 - |O m |/ -JOL + (vFq) 2 ■ This multiplicative depen- 



dence implies that the magnetic field only accounts for regular 
corrections in the form (fisH/Ef) 2 for the q = piece in the 
free energy, i.e., no crossed TH 2 term appears. 



Higher-order terms — As we just found, the second-order 
result for 6x s (T) does not distinguish between short-range and 
long-range interaction, i.e., the specifics of the Coulomb case 
does not show up. There is no guarantee, however, that this 
will remain so beyond the second order. Of particular inter- 
est is whether the divergence of the Coulomb interaction at 
q = affects the expansion of 6x s (T) oc r in powers of r s . 
For this, we computed the corrections to from the third- 
order diagrams. We found that the divergent U(0) still does 
not directly contribute to the spin susceptibility, however the 
prefactor of the linear in T term gets modified due to ver- 
tex corrections to U(2pf). The corrections involve U(2p F ) 
itself and the momentum integrals of the interaction poten- 
tial. The most singular of these corrections accounts for the 
renormalization between U(2p F ) and the spin component of 
the Landau function r s (7r). This renormalization involves 
T 2 W J d 2 pU(k - p)G(p, co)G(p - 2k, to) and yields the mul- 
tiplicative correction to Q in the form 



B = 1 



4 r 
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sin 0/2 + VI -sin6»/2 



sin 6/2 - VT 



■ sin 6/2 
(5) 

where u{8) = {m/2n)U{q = 2p F sin 0/2). One can easily ver- 
ify that for short-range u(8), the integral over 6 converges, but 
for the Coulomb interaction, when u{6) oc 1 / sin(6 ) /2), the in- 
tegral is confined to small 6 and diverges as log 2 . The true 
divergence is indeed cut off by the screening effects. Still, it 
implies that in contrast to a short-range potential, the prefac- 
tor for the linear in T term in the spin susceptibility in the 2D 
Coulomb system is non-analytic in r s . Including screening in 
the usual way, we obtain from (|5jl within logarithmic accuracy 



B = 1 + 



V2r 



[log 2 r s + 0(\ogr s ) + ...], 



(6) 



In a generic Fermi liquid, the full linear in T correction 
to the spin susceptibility is evaluated in the same way as the 
correction to the specific heat |5]. At each order of perturba- 
tion, one selects two bubbles in which one keeps the singu- 
lar frequency dependence of Tl(q, Q m ), and evaluate all other 
bubbles at Qo- The two selected bubbles yield the TH 2 term, 
other bubbles contribute to the prefactor via the renormaliza- 
tion of the 2kp vertex. Extending the analysis in |5| to the 
spin susceptibility we obtain that the prefactor is expressed in 
terms of the spin component of the full quasiparticle scatter- 
ing amplitude at the scattering angle 6 — n. In the explicit 
form, 



m 

— F sp (n) 
m 



(7) 



Here E F = u F p F /2 is the Fermi energy for free fermions, 
m* is the effective mass, and F sp (7r) is the spin component 
of the scattering amplitude. At weak coupling, F sp (n) = 
-mU{2p F )l(2n) - -r s /{2 V2), and Eq. Q reduces to Eq. 
0. In a generic Fermi liquid, F sp (n) - 2„ = o(-l)"(2« + 



3 



1) /sp,n/(l + /sp,«X where / sPi „ are the partial spin components 
of the Landau function. 

If the system is close to a ferromagnetic (Stoner) insta- 
bility, / sp o * -1, and FfJn) can be well approximated by 

4, /(l + /sp.o) 2 Then 6 X s(T) * (T /2E F ) X 2 J X p,uU, where 
Xs = XPau\\(m* lm)l(l + / sp>0 ) is the spin susceptibility in a 
Fermi liquid at T — 0. We did not analyze higher-order terms 
in temperature, but based on the form of 6x s (T), it is tempting 
to assume that 



(8) 



If, on the contrary, the spin component of the scattering am 
plitude is small, but m* jm is arbitrary, as some studies sug 
gest the same consideration yields 

T 

-F 2 M)X v L\v 



x:\t) 



-X~s\T 



0) 



2E f 



(9) 



In silicon inversion layer, £> ~ 6K for typical densities |9]. 
The susceptibility measurements have been reported for T ~ 
2 - AK. xl l (T) measured in units of^p au ii changes by about 
20% between 2K and AK [9]. This would be consistent with 
Eq. however the sign of the measured temperature cor- 
rection is opposite to that in (|SJl. Recent Shubnikov-deHaas 
measurements, however, reported a much weaker, almost un- 
detectable T dependence of X s(T), from which 5 Xs {T) could 
not be extracted 11 411 . This much weaker effect would be 
more consistent with Eq. if we assume that F sp (n) remains 
small. More precise measurements of Xs (T) are clearly called 
for to test our theoretical predictions. 

Extended Landau formalism — We now consider in more 
detail the physics behind the linear in T dependence of the 
spin susceptibility. We argue that this term is actually consis- 
tent with the structure of Landau theory can be viewed as orig- 
inating from the non-analytic structure of the Landau function 
near the Fermi surface. 

We remind that the Landau function is the second varia- 
tional derivative of the energy of the system with respect to 
the distribution function of quasiparticles «o-(p) 



fa* (p, p') 



5 l E 



(10) 



The energy 
field H is | 



Sn tr (p)6n (r (p') 
|ain of a quasiparticle in a weak external magnetic 



6s(p) = -/i B (tH) + Tv f fo* (p, P') ^Wp')-^, 

J os' (2nf 

(11) 

where both 5s and 5n are matrices in the spin space. The g- 
factor of a quasiparticle with momentum p is defined by 



6s(p) = — (crH) , 



(12) 



and the spin susceptibility of a Fermi liquid is expressed as 
the momentum integral over g(p) 

p\ f 8n(p) d d p 



14 r a«(p) <f' P 



(13) 



Eqs. \\ \\ and d!2i determine the integral equation for the 
gr-f actor [7]: 

.(p)=2 + i// 8P (p >P o^(p')^;. d4) 



ds' 



(2n) a 



Here / sp is the spin component of the Landau function: 
v f (P, P') = fJ + trtr'ftp (p, p), where v = m*/n. 

At zero temperature, the integration in d!3> is confined to 
the Fermi surface. At finite T, the quasiparticles are allowed 
to deviate from the Fermi surface. Introducing £ = vp (p - pp) 
and f = vp (p' - pp), and assuming that gr(p) depends on £ but 
not on the direction of p, we re-write Eqs. fl!4i and dl 31 as 

r dn 

=-i-Jl^ (16) 



and 



X 



Here f sp (^, £') is the interaction function averaged over 
the angle <j> between the momenta p and p': / sp (^, £') = 

JT 2 * d<t>K2n). 

At T — 0, only particles at the Fermi surface matter, and 
Eq. Jl 51 yields the well-known result g(T — 0) = 2/(1 + / sp> o) 
For calculations at a finite T, we need the solution at small but 
finite To illustrate the appearance of singular terms, let us 
first study the structure of the finite temperature /-function in 
the limit of weak interactions. The corresponding RPA corre- 
lation energy can be written as follows II 1511 

/de 1 C doj C d 2 Q 
^J^J^ n(W ' q)[y(W ' q) - U(q)] - 

(17) 

In Eq. dl7t . v(q) is the bare Coulomb interaction in two di- 
mensions and V(oj, q) is the dynamically screened interac- 
tion V(«,q) = v(q) [1 - TI(ll), q)v(q)] The polarizability 
YI(cl>, q) is defined as 

n( W ,q) = Ti V f p-. f ^GJs,p)GJs+co,p+q) (18) 
J 2m J (2n) z 

where Ga-(s, p) is the time-ordered Green's function, which 
we write as a functional of the quasiparticle distribution func- 
tion: 



GJs,p) 



MP) 



1 - «ct(P) 



(19) 



s - £o-(p) - iO s - Ea-(p) + iO 

Using the definition fllOi and Eqs. tilt . Jl 8i . and ( 119b . one can 

obtain the quasiparticle interaction function by a straightfor- 
ward evaluation of the derivatives with respect to n^ip) fiH . 
At finite temperatures, we have to allow quasiparticles to de- 
part from the Fermi surface. Keeping this in mind, we find that 
the spin-dependent part of the Landau function is essentially 
the dynamically screened Coulomb interaction, determined by 
the values of momentum and energy transfer of the interacting 
quasiparticles: 



/ S p RPA, (p, P') = -Re V[E(p) - E(p'), p-p']. 



(20) 



4 



In the limit of low temperatures T/E-p <k r s <k 1, we 
can present the spin-dependent part of the Landau func- 
tion /(p, P ') = / sp (£f ;<*) as / sp (^,r;0) = f«g(g,e;4>) + 

/sing <P)< where the "regular" part is due to the statically 
screened Coulomb interaction, and the "singular" part comes 
the singular dynamic screening near 2k p. A small the 
singular part is a small correction to /, es , and can be written as 

/ sing & r ; 0) = ~ [/re g (£ f ; 0)f iWf , f ; 0), (2D 

where n s i ne is the dynamic part of II, given by Eq. at 
oj = E(p)-E(p'). 

For actual calculations of the susceptibility, we will need 
the interaction function averaged over the angle (p. Performing 
the calculations, we find that the singular part is non-analytic 
in the deviations from the Fermi surface. 



/sing Jspfl r — + 

2 l m I tf 



(22) 



where the dots stay for regular terms in £ 2 and (£') 2 - We ex- 
plicitly verified that the non-analyticity in (I22> originates from 
the dynamic 2kf Kohn anomaly. This agrees with our dia- 
grammatic analysis above. 

Substituting the non-analytic part of the Landau function 
into the Eq. Jl 5i for the gr-factor, we obtain the following in- 
tegral equation: 

+ CO 

-E V 

_2 +°° 

This integral equation can be solved by iterations, using the 
zero temperature result g(T — 0) = g* — 2/(1 + / sp ,o) as the 
first approximation. Performing the calculations, we obtain 



9(£,T)= g* 



1+f ^m { *- 2) m 



f 2 _ 



|+21n(l+ e ^ r ) 



(24) 



Using Eqs. dl6> we then obtain for the spin susceptibility 



X = 



9j ^°m 8E F 



(25) 



This agrees with the result of the diagrammatic treatment, 
Eq. 0, for the case when the full scattering amplitude F sp (n) 
can be approximated by its zeroth partial component F^„ « 
fspo/(^ + / s P.o) 2 - This approximation is implicit in the RPA 
formalism. In a more generic analysis, one indeed should re- 
cover the full scattering amplitude. 

We see therefore that the linear in T correction to the spin 
susceptibility can be understood as originating from the non- 
analytic momentum dependence of the Landau function at 



small deviations from the Fermi surface. This non-analytic 
momentum dependence / s ; ng oc |£ — £'| is the fundamental 
consequence of the dynamic 2pp Kohn anomaly in a generic 
Fermi liquid. 

To conclude, we considered the temperature dependence of 
the spin susceptibility in a 2D electron system with Coulomb 
interaction. We found that the leading temperature correc- 
tion 3xs(T) is linear in T and comes from the 2pp singular- 
ity in the dynamical response function (the dynamic Kohn 
anomaly). The origin of the effect is the same as in systems 
with short-range interaction. However, for Coulomb inter- 
action, the prefactor of the 0(7") term is itself non-analytic 
in r s . We also analyzed the emergence of the 0(T) term in 
5xs(T) by extending the Landau formalism to finite T. We 
demonstrated that within this approach, the non-analytic tem- 
perature dependence of the susceptibility originates from the 
non-analytic momentum dependence of the Landau function 
at small deviations from the Fermi surface. The experimental 
verification of the linear in T dependence of the spin suscep- 
tibility is clearly called for. 
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